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ABSTRACT 

We present 2" - 10" imaging of eleven transitions from nine molecular species across the nuclear 
bar in Maffei 2. The data were obtained with the BIMA and OVRO interferometers. The ten 
detected transitions are compared with existing CO isotopologues, HCN, CS and millimeter continuum 
data. Dramatic spatial variations among the mapped species are observed across the nuclear bar. 
A principle component analysis is performed to characterize correlations between the transitions, 
star formation and molecular column density. The analysis reveals that HCN, HNC, HCO+ and 3 
mm continuum arc tightly correlated, indicating a direct connection to massive star formation. We 
find two main morphologically distinct chemical groups, CH3OH, SiO and HNCO comprising the 
grain chemistry molecules, versus HCN, HNC, HCO+ and C2H, molecules strong in the presence of 
star formation. The grain chemistry molecules, HNCO, CH3OH and SiO, trace hydrodynamical bar 
shocks. The near constancy of the HNCO/CH3OH, SiO/CHgOH and SiO/HNCO ratios argue that 
shock properties are uniform across the nucleus. HCN/HCO+, HCN/HNC, HCN/CS and HCN/CO 
ratios are explained primarily by variations in density. High HC0"'"/N2H+ ratios are correlated with 
the C2H line, suggesting that this ratio may be a powerful new dense photon-dominated region (PDR) 
probe in external galaxies. C2H reveals a molecular outflow along the minor axis. The morphology 
and kinematics of the outflow are consistent with an outflow age of 6-7 Myrs. 

Subject headings: galaxies: individual(Maffei 2) — galaxies: starburst — galaxies: ISM — radio lines: 
galaxies — astrochemistry 



1. INTRODUCTION 

Molecular abundances and excitation are sensitive 
probes of p hysical processes taking plac e in molecular 
clouds (e.g.. Ivan Dishoeck fc Blakd ll998^. Chemistry is 
pivotal both in controlling gas cooling and ionization 
balance as well as in revealing changing cloud physical 
conditions resulting from shocks and feedback from mas- 
sive star formation and galactic structure. Millimeter 
telescopes are now capable of surveying many molecular 
transitions in other galaxies, opening up chemical studies 
to the rest of the Universe. 

Single dish millimeter line surveys show that spec- 
tra of external galaxies are extremely rich in molec 
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Costagliol a et al 
2011bll . For nearby galaxies, giant molecular cloud 
(CMC) scale chemistry is within reach of current mil- 
limeter interferometers. In the fir st pape r of this series 
(Meier & Turner 2005; hereafter iMTOSf ). we exploited 
this richness to carry out the first chemical imaging sur- 
vey of the nucleus of an external galaxy with spatial reso- 
lution high enough to separate individual GMCs. Images 
of eight astrochemically important species were obtained 
in the nearby spiral IC 342, revealing strong chemical 
differentiation among the GMCs. A principal compo- 
nent analysis demonstrated that the chemical variation 
is strongly correlated with galactic structure such as bar 
and arm features, as well as the presence and evolu- 
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tiona ry state of star forming regions (jMTOa lM eier et aU 

Eoni). 

In this paper we extend the analysis to the nucleus of 
a second nearby spiral galaxy, Maffei 2. We have two 
main goals for this survey. The first is t o test the gen- 
eralizability of the results discovered in IMT051 Maffei 
2 and IC 342 are at the same distance and their nuclear 
morphologies are similar. Both have two molecular spiral 
arms terminating on a central ring, with the central ring 
being the site of intense, young, massive star formation. 
It is reasonable to expect that the two nuclei should ex- 
hibit similar chemistries if the results of lMTOSi for IC 342 
are applicable to barred spiral nuclei in general. 

The second goal is to investigate the chemistry of Maf- 
fei 2 with partic ular attention focused on its unique 
features. IMT05I find that in IC 342, shocks appear 
to be an important influence on the nuclear chemistry 
and speculate that the shocks are associated with or- 
bital reso nances in a bar potentia l. The bar in IC 342 
is weak (jCrosthwaite et all |2001| ): however, and there 
remain other explanati ons for the nuclear morphology 
(jSchinnerer et al.) I2008D . It is important to follow up 
with an object that has a well established, strong bar. 
Kinematic modeling leaves little doubt that the nucleus 
of Maffei 2 hosts a strong, highly defined nuclear bar 
of 320 PC radius (Meier et al. 2008, hereafter iMTHOSi : 
iKuno et "aLll2008| ) . Because of the strength of the bar, we 
expect that shock chemistry should be more important 
than in IC 342. Maffei 2 also has a some what more in- 
tense and extended starburst than IC 342 (jTurner fc Hoi 
1983, 1994), so photon-dominated region (PDR) effects 
may be more prominent than they are in IC 342. 

2. OBSERVATIONS 
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To pursue these goals, we obtain aperture synthe- 
sis observations of eleven important molecular transi- 
tions over the central R~0.5 kpc of the nuclear bar 
in MafFei 2, for comparison with CO isot opologues 
(|Hurt fc T urner T99T: 'MTH08), HCN(l-O) (IMTH08I) 
and CS p^uno ct al. 2008) . Table [T] lists molecular pa- 
rameters for each transition. With the Owens Valley 
Radio Observatory (OVRO) millimeter interferometer 
we observed three lines, HCO+(1-0), HNCO(5o5-4o4) 
and SiO(l-0; v =1). HCO+(1-0) and SiO(l-0; v 
1) were observed between 1999 January 28 and 199 9 
March 29 along with the HCN(l-O) dataset (|MTH08l ). 
HNCO(5o5-4q4) was observed along with "CO(l-O) and 
0^^0(1-0) between 1998 October 19 and 1999 January 
5 (Table [1]). The OVRO interferometer consisted of 
six 10.4 meter antennas equipped w ith SIS receivers 
(jPadin et al.lll991l:lScoville et al.lll99"l . Transitions were 
observed in L, H and UH array configurations except for 
the HNCO, which was observed in L and H, providing 
spatial resolutions of 2. 5 (Table [1]). Data were cal- 
ibrated using the in-house OVRO calibration package, 
MMA. Phase calibration was done by observing the point 
source 0224-1-671 every 20 minutes. Absolute flux cali- 
bration was based on observations of Neptune or Uranus 
with 3C 273, 3C 84 and 3C 454.3 as supplementary flux 
calibrators. Based on the derived fluxes and flux histories 
of these secondary flux calibrators we estimate that the 
absolute fluxes are good to 10 - 15%. Two pointings were 
mosaiced to cover the nuclear bar and therefore maps are 
corrected from the primary beam response. 

The eight remaining transitions were observed with 
the ten element Berk eley - Illinois - Maryland Associ- 
ation (BIMAj3 array ()Welch et all [1996) between 2002 
September 08 and 2003 March 18 (Tabled]). The eight 
spectral lines were observed in two sets of spectrometer 
configurations. C2H(N = 1 - 0, J = C2H(N = 1 - 

0, J = i-i), SiO(2~l;w = 0), HNC(l-O) and HC3N(10- 
9) were observed together, as were CH30II(2fe — 1^), 
C34S(2-1) and N2H+(l-0). Data were obtained in B, C 
and D configurations for all transitions giving beamsizes 
of ^ 7 — 10 depending on frequency and sensitivity. 
Calibration and data reduction was done in MIRIAD. 
Phase calibration was also done by observing the same 
point source 0224-^671 (0228+673) interlaced with the 
source every 24 minutes. Absolute flux calibration was 
done by observing W 3(011). Based on the derived fluxes 
and flux histories of 0224-1-671 we estimate that the un- 
certainties in absolute flux are 10 - 15%, similar to the 
OVRO data. The - 130" FWHM primary beam of the 
BIMA antenna is large enough such that only one point- 
ing was required. No corrections for primary beam at- 
tenuation have been applied. 

Image analysis was done with the NRAO AIPS pack- 
age. In making the integrated intensity maps emission 
greater than 1.2cr was included. Continuum emission has 
not been subtracted from the maps since the 3 mm con- 
tinuum peak is below la. 

3. RESULTS 

Operated by the University of California, Berkeley, the Univer- 
sity of Illinois, and the University of Maryland with support from 
the National Science Foundation. 



3.1. A Sketch of the Nuclear Region of Maffei 2 

The Sb galaxy MafFei 2 is one of the nearest and bright- 
est extragalactic molecular li ne sources, at a dist ance 
of only 3.3 Mpc (1" ~15 pc, iFingerhut et~al] [20071 see 
discussion in MTH08). A strongly barred and some- 
what asymmetric galaxy, Maffei 2 has undergone a re- 
cent interaction with a small companion that has driven 
a large quantity of gas into the nucleus, building up 
a compact central bulge seen in the NIR ([Hurt et al.l 
119931 119960 . This bar-like compact bulge, has further 
driven the molecular gas into a central molecular zone 
of radius R ^ 300 pc that is barred and h ighly inclined 
(jlshiguro et al.ll 19891: [Hurt fc Turnerill99lD . This nuc lear 
bar, containing M(H2)- 2 x 10^ Mq of gas (IMTH08I ). is 
structurally remini scent of a miniature version of large 
scale gaseous bars (| Athanassoulal [l992[ ) (Fig. [T]). Using 
2 cm and 3 mm radio / mm continuum to trace ion- 
izing radiation it is seen that the intersections of bar 
arms and central ring are sites of intense localized star 
formation (SFR ^ 0.04 Mq yr"^ per CMC or a total 
of Lqb ~ 10* L^: iTurner fcHol[T99l [Tsai et al.l[2?}06l : 

IMTH08D . Averaged over the inner ~ 20 radius, gas 
consumption timescales decrease below 100 Myrs, mak- 
ing Maffci 2's nucleus a starburst. 

A central molecular ring of radius ^ 7" (110 pc) is 
clearly resolved in gas kinematics, although t he ri ng is 
less apparent in integrated intensity images ([MTH08). 
GMCs D, E, and F (cloud designations following MTH08 
are shown in Figs.[T|) are located along the eastern side of 
the nuclear ring. Bar arms (including GMCs B and C in 
the north and G and H in the south) extend off the lead- 
ing edges of the ring. The bright GMCs of the central 
ring, typically ~60-80 pc in extent, and with masses of a 
few 10® Mq, delineate the sites where gas flowing in along 
the bar arms piles up at the inner ring. Dense gas traced 
by HCN preferentially picks out the arm-ring intersec- 
tions. Gas not incorporated into the dense component 
at the starburst locations is tidally shea red into a mod- 
erate density, smoothly distributed ring (|MTII08l ) . Mas- 
sive star formation is largely absent in the western side of 
the ring (MTII08). Based on HCN, the fraction of dense 
gas is significantly lower towards the bar ends. Densities 
derived from CO tend to be fairly constant {nH2 ~ lO'^) 
over the inner bar, with the clouds associated with young 
star formation slightly warmer {Tk 30 - 40 K) than 
the others (T^ ~ 20 K). Gas clouds on the quiescent 
western side of the nuclear ring are slightly cooler yet 
(|MTH08. ). 

3.2. Molecular Abundances in Maffei 2 

Column densities are determined assuming optically 
thin emission, and LTE: 

/ SkQeW. N 

where S^i, g and are the line strength, degeneracy 
and upper energy of each state, respectively, and Tex 
is the excitation temperature associated with the tran- 
sition. Given that we have mapped only one transition 
of each species, corrections for background radiation and 
opacity have been ignored. Column densities are sensi- 
tive to Tex through the partition function, Q, and the 
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energy of the upper state. The asymmetric tops (HNCO 
and CH3OH) are more sensitive to temperature changes 
than are the hnear rotors. Changes in gas density also af- 
fect excitation, particularly for molecules with high crit- 
ical densities (HNC, N2H+ and HC3N). Densities have 
been at least p artially constrained by observation in Maf- 
fei 2 ()MTH08f ). We adopt a T^x = 10 K for the molecular 
transitions here. This is similar to what is derived from 
the "CO d ata (IMTHOST l and is the same as assumed 
for IC 342 (|MT05l ). permitting easy comparison. How- 
ever Tex = 10 K may slightly underestimate excitation 
at the massive star forming regions. Measured line inten- 
sities, peak temperatures, line centroids and line widths 
for the GMCs are reported in Table [T] Fractional abun- 
dances (X(mol) = Nmoi/NHa) Calculated from the inten- 
sities given the above assumptions are listed in Table [U 
based on molecular parameters in Table [TJ 

Fractional abundances also require an H2 column den- 
sity, N(H2). N(H2) is most easily obtained from the 
CO (1-0) brightness and an empirical Galactic conver- 
sion factor, Xco- However, Xco overpredicts N(H2) in 
nearby galaxy centers, including Maffei 2, by factors of a 
few (jMTHOS). An alternative measure of N(H2) can be 
obtained from the CO isotopologues, which tend to be 
optically thin. We use the ^^CO(l-O) integrated inten- 
sity and [H2/13CO] ^ 7.0 X 10^ = 60 and 
[CO/H2] = 8 .5 X 10~^ ) when calculating the H2 column 
densities (see IMTH08I . for a detailed discussion of N(H2) 
and its uncertainties in Maffei 2). 

Fractional abundances for the species below are esti- 
mated to be uncertain to at least a factor of three, al- 
though the relative column densities — that is, the rela- 
tive spatial distributions within the nucle us of th e dif- 
ferent molecules — should be more reliable. IMT05I has a 
detailed discussion of the uncertainties associated with 
the abundance estimates. Table [1] records for reference, 
the change in derived column densities if Tex is changed 
from 10 K to 50 K, more typical of the large beam gas 
kinetic temperatures. 

3.3. The Molecules of Maffei 2 

We now discuss the morphologies of each species. Spec- 
tra taken at selected locations across the nucleus with 
8" apertures are shown for each transition in Fig. [2l 
Fig. [3] displays maps of the transitions detected with 
OVRO and Fig. H displays those detected with BIMA. 
In the subsequent analysis we also includ e 3mm co ntin- 
uum emission and CO isotopologues from|MTHO|. The 
millimeter continuum in Maffei 2 is do minated b y free- 
free emission rather than dust emission (|MTH08n . 
SiO — Silicon Monoxide: The J = 2-1 rotational line 
of both the v — Q and 1 vibrational states were ob- 
served. The V = 1, J = 2 - 1 tran sition is a maser 
in the Galaxy (jSnvder fc Buhll I1975D . but we do not 
detect it in Maffei 2. The thermal v = Q transition 
of SiO is tentatively detected toward GMCs B, F and 
G with peak temperatures over the 8" aperture of 35 
mK (Fig. S]). The integrated intensity map is noisy 
but consistent in morphology to HNCO and CH3OH, 
extended along the southern arm. Implied SiO abun- 
dances reach ~ 1 x 10^^". This abundance is lower than 
estimated for IC 342 and the nearby starburst nuclei, 
NGC 253 and M 82 (Garcfa-BuriUo ct al.. .2000.. .20011 



lUsero et al.ll2006 D. but t ypical of diff use Galactic clouds 
(|Greaves et Wr imQ: T mne^ll998bD away from molec- 
ular outflows (Martm-Pintado et al.l I1992D ; it is much 
larger than what i s found for Galactic dark clouds (e.g., 
iZiurvs et al1[T98l . 

C2H — Ethynyl: We detect both the J = f - 5 and | - | 
fine structure components of the N = 1-0 line. C2H is 
brightest towards the northern star forming ring (CMC 
D+E), peaking at 0.11 K, and weakens significantly along 
the bar arms (Fig. H]). Towards GMCs F and G, C2H 
avoids the column density peaks traced in ^^CO(l-O). 
Line widths are significantly broader than those seen in 
the other weaker lines, so the hyperfine structure of each 
fine structure component is being marginally resolved 
here. Evidence for the most blueshifted N = 1 - 0, J 
— i — i,F = l- hyperfine component is seen towards 
the brightest GMCs (see Fig. [5]). There is a plume of 
emission along the minor axis in the | — ^ line. From 
the relative intensities of the fine structure lines we can 
determine the C2H opacity. The ratio of the bright- 
ness temperature of the two fine structure transitions, 
-^(i ~ ~ h) ranges from 1.2 - ^2.8, with lowest val- 

ues toward the regions with the brighter emission. For 
optically thin, LTE excitation a this ratio should be 2.0. 
Therefore C2H opacities are significant along the cen- 
tral ring, with implied values for the § — 5 transition 
of TccH — 1-5. Along the arms tcch <1 except at 
regions around CMC G away from the ^'^ CO (1-0) col- 
umn density peaks where tcch 3. Hence, despite its 
lower signal-to-noise, the ^ — 5 transition probably rep- 
resents a closer picture of the true C2H abundance dis- 
tribution. Fractional abundances assuming C2H is op- 
tically thin peak at ~ 2.5 x 10^^. If we account for 
opacity, the abundance rises to ^ 4.8 x 10^^. These are 
similar to those seen towards the FDR site in IC 342 
(|MT05[ ) and at the high end of the range seen in Galac- 
tic co res (Wootten et al. 1980; Huggins et al. 1984 Hmt^ 
Il983t [Turn er et al. 1999; L ucas fc Lis zt 2000) . 
HCO+ —Formyl Ion: HCO+(1-0) is the brightest line 
of this sample, and is imaged at high resolution (^2 ). 
The emission is dominated by GMCs D-t-E and F, with 
line peaks of 2.0 K ( 0.58 K o ver the 8 aperture) (Fig. 
[31). Like HCN(l-O) (|MTH08D . HCO+(1-0) is primarily 
confined to the central ring region. There is a second 
blueshifted velocity component at ~-135 km s^^ (LSR) 
along the northern portion of the bar. HCO+(1-0) abun- 
dances in Maffei 2 are between 2 - 8xl0~^, but these 
(along with those of HCN and HNC) should be consid- 
ered lower limits since these very bright transitions may 
have significant optical depths (see section l4.1.ip . These 
agree well with the HCO"*" abundances i n most Galac- 
tic d e nse clouds (eg . iPirogov et al.l 119951 : iLucas fc Lisztl 
119961 : lTurneiill995aD and single-dish e xtragalactic mea- 
surements l|Nguven-0-Rieu et al.iri992D . 
HNC — Hydrogen Isocyanide: The J = 1-0 line of the lin- 
ear molecule HNC is bright over the whole nuclear bar, 
with peak temperatures '^0.25 K (Fig. [4|). At this reso- 
lution, the morphology matches closely both ^'^CO(l-O) 
and HCN(l-O), being dominated by the bright clouds, D 
+ E and F, but HNC emission is seen from every cloud in- 
cluding the non-nuclear bar CMC, A. Towards the north- 
ern side of the bar HNC emission is especially bright at 
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higher blueshifted velocities, similar to the component 
seen in HCO+(1-0). HNC a bundan ces are 1 - 2 x lO'^, 
in agreement with IC 342 ()MT05[ ). These abundances 
are intermediate between typical massive star forming 
cores and Galact i c dark cloud s (Woottcn ct alJ |1978l: 
Blake et alj 119871: IHirota et ahl J998.; N umnielin et al l 
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HC3N — Cyanoacetylene: We mapped the J = 10-9 rota- 
tional line of HC3N. This molecule has the largest electric 
dipole moment and the highest upper energy state of the 
(thermal) sample (Table [1]). It is thus not surprising that 
we observe HC3N to be faint and spatially confined (Fig. 
2]). To improve S/N, the HC3N map was smoothed to a 
slightly larger beam size, but even then peaks at only ~ 
50 mK towards two unresolved locations at the intersec- 
tion of the ring and bar arms (C and F). Abundances are 
~ iq -Q, similar to those in IC 342 (MTO|; iMeier et alJ 
I2OIID , somewhat higher than in cold Galactic clouds but 
much lower than the localized Galactic hot core values 
(> 10"'^; eg. lMorris et al lll976HVanden Bout et"all ll983l: 
iChung et alM ggi: dc Vicente et al.ll2000t ). 
N2H'^ — Diazenylium: In the Galaxy, N2H+ is regarded 
as a tracer of quiescent gas. We observed the J — 1-0 ro- 
tational state of this linear ion, which has hyperfine split- 
ting much smaller than the observed linewidths. As is the 
case in IC 342, N2H+ is bright and extended across the 
nucleus (Fig. |4]). Emission peaks towards CMC F with 
a peak brightness temperature reaching ^0.15 K. GMCs 
D and E, the closest clouds to the northern star for- 
mation region, show weak emission in N2H+: X(N2lI+) 
reaches ^ 7 x 10^^° at CMC F and drops by nearly 
a factor of three towards the 3 mm continuum peak. 
The highest abundances are similar to what is seen in 
IC 342 and de nse, q uiescent cores in the Galactic disk 
(jWomack et al.lll992l : iBenson et al.lll998D . 
(?^S —Carbon Monosulfide: This is an isotopologue 
of the dense gas tracer CS. The J = 2-1 rotational 
transition of C^'^S is tentatively detected only towards 
GMCs A, C and F, at < 30mK (Fig. H). The 
tentative detections imply A"(C'^'*S) abundances of at 
most a few xlO^^", typical of Galactic dense cores 
and diffuse/translucent clouds (for a C'^^S/C'^'^S isotopic 
abundance of 23; eg. 'Nvman 1984; Drdla_ct al. 1989; 
Wang et al.|[T9" 93: Chin ct al. 1996; Mor ata et all Il997t 
Lapinov et alTri998.: .Lucas fc Liszt.. 2002. ). In the GMCs 



near the northern star forming region abundances are 
an order of magnitude lower. Evidently C'^'*S(2-1) is 
not strongly enhanced anywhere across the nucleus, and 
especially not at the massive star formation sites. Lim- 
its on the CS/C^"'S isotopologue ratio assuming optically 
thin CS ranges from 4 to >14 based on the CS(2-1) line of 
IKuno et all (pOOl . The tentative detections of CS/C^^S 
~ 4 suggests opacities significantly greater than unity in 
CS(2-1). 

CH3OH — Methanol: We observed the low energy, 
blended set of 2i - liE, 2q - IqE, 2o - loA-f and 
2_i — l_iE thermal transitions of CII3OII (hereafter la- 
beled the 2fc — Ifc transition). We cannot resolve these 
transitions in this group, given the source linewidths. 
The 2k — Ik transition of CII3OH is very bright across 
the nucleus, with peak temperatures reaching 0.28 K 
(Fig. |4]). Of the species observed in this paper only 
IICO+(1-0) is brighter. Remarkably though, the mor- 



phology of CII3OH is completely different from what is 
seen in most of the other lines, including ^"^00(1-0) and 
IICN(l-O). CH3OII peaks strongly towards the southern 
nuclear bar end (G and H) and the southern nuclear bar 
arm / ring intersection region (F). By contrast, emission 
is barely detected near CMC D-l-E. The emission favors 
the right side of the southern arm, which is the trailing, 
upstream side. At the bar ends methanol abundances 
peak at ~ 1.5 x 10~^. Towards GMC D abundances drop 
by a factor of 5. The peak arm abundances are a factor 
of two higher than those of IC 342 (MT05). Galactic 
methanol abundances can re ach these levels in srnall lo 
calized hot cores regions (eg. iKalenskii fc Sobolevl 11994 ; 
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Turned 1199831 (Nummelin et al.l l2000t iMinier fc BoothI 
20021) . CH3OH velocities appear consistently blueshifted 



with respect to most of the other transitions. In the 
north the velocity shift is ~ -10 km s~^, growing to ~ 
-40 km along the southern bar end. The shift to- 
wards the northern cloud is not significant, however the 
shift seen towards GMCs F and G appears real. This 
may result from an increased population in the higher 
excitation blueshifted lines that make up the quadruplet 
or somewhat different kinematics due to the presence of 
shocks (see section H.4.ip . 

HNCO — Isocyanic Acid: We observed the K_i — tran- 
sition of the J = 5-4 rotational state (5o5 — 4o4) of this 
prolate, slightly asymmetric top. In the higher resolution 
OVRO maps, HNCO peaks at GMCs F and G, with faint 
emission along the bar connecting the two clouds (Fig. 
[3]) . When smoothed to match the resolution of the BIMA 
data, HNCO has a very similar morphology to CH3OH. 
Over the 8 aperture used for the spectra in Fig. [21 
HNCO has a peak brightness of 0.25 K, whereas over the 
40 pc resolution achieved in the OVRO data peak bright- 
nesses reach 1.0 K. No emission is detected towards GMC 
D in the high resolution map. HNCO abundances reach 
X(HNCO) ~ 9 X 10"^ along the southern nuclear bar. 
These abundances are typical of those observed on one to 
two orders of ma gnitude smaller scales in Galactic mas- 
sive dense cores i)Zinchenko et al.ll2000[ ) and are even a 
factor of 3 larger t han the HNCO enhanced shock regions 
in IC 342 (lMT05r ). CH3OH and HNCO in Maffei 2 have 
some of the most pronounced morphological differences 
from CO yet seen in extragalactic sources. These two 
species illustrate the importance of high spatial resolu- 
tion interferometric observations of chemical species. 
Other Transitions: Table [1] presents the limits (2cr) for 
other potentially detectable transitions within the ob- 
served bandwidth. These include the ^'^C substituted iso- 
topologues of HNC(l-O) and HC3N(10-9) and two tran- 
sitions of the radical HCO. None of these transitions are 
detected. In the case of the "'^^C isotopologues only the 
HNi3c(i-0) limits are at all constraining (Table[l]). To- 
wards the central ring and southern bar end la HNC(1- 
0)/HNi3c(2-0) line ratio lower limits are >10-16, which 
in these locations are consistent with or slightly a bove the 
corresponding CO isotopic line ratios (jMTH08D . There- 
fore the HNC(l-O) opacities are lower than ^^00(1-0) 
opacities. 

To facilitate the exposition, we omit the transitions 
for each species where possible in the remainder of the 
text. We discuss transition dependent effects specifically 
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where important. 

3.4. Understanding the Molecular Maps: Principal 
Component Analysis 

The maps give a picture of the molecular cloud astro- 
chemistry in the nuclear region of Maffei 2. To inter- 
pret these maps, we need to establish similarities and 
differences between the molecules. The correlations re- 
veal trends in the chemistry, which in turn, can suggest 
the physical driving forces responsible for these charac- 
teristics and how they change across the nucleus. 

As in MT05, we quantify morphological corre- 
lations amongst the images via a principal com- 
ponent an alysis (PCA). PGA is a common tech- 
nique l^eg. iMurtaee fc Heckl Il987t iKrzanowskil 119881 : 
lEveritt fc Dunnl 120011 : IWall fc ,Tenkins ll2003D used to re- 
duce the dimensionality of a dataset. The PCA sim- 
plifies the picture of molecular distribution, reducing a 
large amount of information to a few images, providing 
an excellent framework within which to study the com- 
plex variations in molecular properties in Maffei 2. For 
a general descriptio n of PCA applied to mult i-transition 
molecular maps, see IHeyer & SchloerB (|T99l or for spe- 
cificapplications to extragalactic molecular maps, see 
IMT051 The PCA is dependent on the particular set of 
data given to it; here, we are analyzing the correlations 
of 3 mm lines of heavy molecules, so the results we obtain 
pertain to the cool and relatively dense cloud population 
in Maffei 2. 

To calculate the principal components for Maffei 2, 
the line maps were remade to the same geometry and 
beamsize (that of HNC; Table [1]), normalized and mean- 
centered. Thirteen maps, includ ing CO (1-0), HCN(1- 
0) and 3 mm continuum from IMTH081 plus the ten 
transitions of the nine species discussed in section 13. 3[ 
were sampled at 1 intervals over the central 38 x 56 , 
making up a 27664 element data-space. The algorithm 
used to calculate the principal components is that of 
IMurtage fc HecS ([1987). 

The results of the PCA are tabulated in Tables [T] and 
[TJ Table [T] shows the matrix of pair- wise correlation coef- 
ficients between all of the thirteen maps. Table [T] records 
the projection of each map onto the seven modeled PCs. 
Fig. [5] displays the maps of the first three (most im- 
portant) principal components (PCs). These three PCs 
account for 80 % of the variation present in the data set. 
The projections of each map onto the first three PCs are 
shown vectorially in Fig. [51 

Inspection of Figure [5] reveals that PC 1 is most 
strongly correlated with all maps. It represents a mix 
of the ^^CO and HCN maps, implying that PC 1 is a 
good density-weighted column density map of the low 
excitation molecules in Maffei 2. This was also found 
in IC 342; however, here ^^CO projects along PC 1 more 
closely than HCN, hence PC 1 is nearer to a pure column 
density map than its analogue in IC 342. All maps ex- 
cept the low S/N ratio C^'^S show significant correlation 
with PCI. 

PC 2 characterizes the dominant variation seen once 
the distribution of column density has been accounted 
for. PC 2 distinguishes two major distinct groups of 
molecules, those peaking towards the norther star form- 
ing region and those peaking towards the southern nu- 



clear bar end. Species with large positive projections 
on PC 2 (star formation peakers) include HCO"*", HCN, 
HNC, C2H and of course 3 mm continuum, which is free- 
free emission. Species with large negative projections 
(southern bar peakers) include HNCO, CH3OH and to 
a lesser extent N2H+ and SiO. PC 3 has a quadrupo- 
lar morphology dominated by SiO and C^^'S. PC 4 and 
beyond are at best marginally significant so we do not 
discuss them in detail. 

The correlation matrix resulting from the PCA (Table 
[1]) shows in more detail the correlations between various 
species. The following trends are noted. 

(1) HCN, HNC, HCO~^ and 3 mm continuum are all 
extremely tightly correlated, indicating a close associa- 
tion with massive star formation. At the resolution of 
the PCA HCO"^ and HCN are effectively perfectly cor- 
related. At three times higher spatial resolution sub- 
tle differences become apparent (Fig. [7|). HCN is the 
map most tightly correlated with the 3 mm continuum. 
Therefore, the (massive) star formation rate versus HCN 
line luminosity relatio n seen in single-dish surveys (e.g. 
iGao fc SolomonI [20041 ) per sists to a t least 100 pc scales 
in Maffei 2, as noted in (lMTH08f ). However, the cor- 
relation between HC3N(10-9) and 3 mm continuum is 
weaker than found for IC 342. 

(2) CH3OH and HNCO are very well correlated with 
each other and are strongly anti- correlated with the mas- 
sive star formation. This leaves little doubt that on these 
CMC spatial scales the chemistries of these two species 
are either directly coupled or are mutually coupled to 
a specific chemical mechanism. Moreover these species 
show little correlation with 3 mm continuum and hence 
star formation is not relevant to establishing the mor- 
phologies of these species. Nor are these species corre- 
lated with the standard dense gas tracers, HCN, HCO+ 
and HC3N. Fig. [8]displays the position-velocity diagrams 
of "CO, C2H and HNCO. The strong differences be- 
tween the starburst tracers and CH3OH or HNCO carry 
over to the kinematics. The 'parallelogram' originating 
from the central ring is nearly absent in CH3OH and 
HNCO but remains pronounced in C2H. 

(3) SiO and especially C^'^S are less tightly correlated 
with the other species. This may in part be due to their 
low S/N ratio and that some of the structure that ap- 
pears in these maps are nois e. However C^'^S was also 
anomalous in IC 342 ()MT05D . 

4. DISCUSSION 

The pronounced variations in spatial morphology ob- 
served for the different molecules across the central half 
kiloparsec of Maffei 2 are caused by either differences 
in emissivity (excitation) or abundances (chemistry). In 
this section we discuss the possible influences on the 
molecular emission that could lead to the observed vari- 
ations. 

4.1. The Dense Gas Component: Diagnostic Ratios 

The CO /HCN ratio is widely used to study the rela- 
tive fraction of dense gas in galaxies. IMTHO^ used the 
CO/HCN ratio to identify the sites of high dense gas 
fraction and concluded that the xi — X2 orbital intersec- 
tions and the eastern (star forming) portion of the central 
ring are the locations of the high dense gas fraction. The 
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fraction of dense gas drops towards the western ring and 
moving out along the bar arms. 

To study the nature of the dense component directly, 
line ratios between dense gas tracers are required, 
since CO does not trace dense gas. Line ratios be- 
tween bright, high critical density species, such as 
HCN, HCO+, HNC, HC3N and CS are useful di- 
agnostics of the physical and chemical properties of 
the dense component of the ISM. Much observational 
and theoretical effort has been dedicated to under- 
standing these diag nostic species in nearby starbursts 
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difficult since most of these studies have been done at 
low spatial resolution, either with single-dish telescopes 
or in distant sources. These ratios mix many interstellar 
environments into one beam. 

Now that we have access to an inventory of dense 
gas tracers at 30-100 pc resolution, we investigate 
changes in the HCN/HCO+, HCN/HNC, HCO+/N2H+, 
HNC/HCO+ and HCN/CS line ratios on GMC scales 
across the nucleus of Maffei 2. Fig. [7] displays four of 
these line ratios. Line ratios are tabulated in Table [TJ 
Organized variations in all ratios are observed on these 
scales and are discussed in terms of the CO/HCN ratio. 

4.1.1. HCN/HCO'^: Tracer of Intermediate Gas Densities 

The HCN/HCO+ ratio relates two of the brightest lines 
that trace dense gas. It has become a popular diagnos- 
tic of gas with densities intermediate between CO and 
HCN since HCO+ has a critical density nearly an or- 
der o f magnitude lower than H C N fe.g.. iKohno et al.l 
20011: iGracia-Carpio etHI l2006t llmanishi et al.l 120071 : 



Krins et al.ri2008f ) . HCO+ also has a formation route 



that maintains significant abundances in relatively dif- 
fuse gas (see below) , and it is an ion making it more sensi- 
tive to collisions with electrons fe.g. lPapadopoulosll2007D 
(see Table [1]) . As a result HCO+ can reinain a bundant 
and excited in diffuse gas (jLucas fc LisztlfToOfih . Theo- 
retical modeling finds that HCN/HCO"'" is a strong func- 
tion of density for a given column density in the presence 
of ionizing radiatio n, with ratios of unity for densities 
of^lO^ cm-3 fe.g. iMeiierink et al.|[2007t lYamada et al.l 
[2003). 

There is, however, the complication of opacity. Despite 
its lower critical density HCO+ has a higher Einstein 
Ay , implying higher radiative efficiencies than HCN (and 
HNC). For a given abundance per line width, HCO"*" ex- 
hibits a line opacity nearly twice that of HCN and HNC. 
For extragalactic studies it is usually assumed that HCN 
abundances are significantly higher than HCO+. How- 
ever, for the LTE abundances calculated in § 13.21 and 
IMTH081 HCO+ would have approximately twice the ob- 
served abundance of HCN, contrary to what is commonly 
assumed. These abundances have large uncertainties but 
the two lines have similar energy levels. The similarity in 
excitation requirements of the two lines makes it difficult 
to reverse this situation purely by changing gas temper- 



atures. If the LTE abundances are correct then HCO+ 
line opacities must be at least as large as the presumably 
large HCN line opacities. If opacities in both lines are 
greater than unity the diagnostic power of HCN/HCO+ 
is limited. 

Towards the nuclear bar in Maffei 2 the overall mor- 
phologies of HCO+ and HCN are similar. However, 
detailed differences in the two transitions are revealed 
in their ratio. HCN/HCO^ ranges from 0.6-2.2 over 
30 pc scales in Maffei 2. These values cover much of 
the rang e sampled in the s i ngle -dish su rveys of many 
galaxies (IKohno et al.l 120011 : iGra cia-Car pio et al.l 120061 : 
iKrips et al.ll2008[ ). Both HCN pTHOSi ) and HCO+ fall 
off more rapidly with galactocentric distance than does 
^^CO(l-O). HCO+ decreases more slowly than HCN in 
the north (GMC B) but more rapidly along the southern 
bar arm. HCN is elevated relative to HCO+ on the active 
eastern side of the central ring, GMC G and a localized 
cloud just south of GMC F, while HCO+ is enhanced, 
relatively speaking, in the quiescent, western ring region. 
We detect clouds with CO/HCN vs. HCN/HCO+ values 
that reside in the "AGN dominated" portion of param- 
eter space set by global single-dish studies where clearly 
there is no AGN (jKohno et al.ll200"TI : IGracia-Carpio et al.l 
l2006f ) . Similar statements apply for " star formation dom- 
inated" parameters. As a result interpretation of single- 
dish values must be done with great care. 

The enhancement of HCO"*" towards the western ring 
and HCN towards the northeastern ring makes sense 
in the conventional critical density context, since one 
expects the denser gas to be more closely associated 
with the massive star formation. However, the elevated 
HCN/HCO"'" ratio towards the southern bar end is un- 
usual. A decreasing HCN/CO ratio suggests the fraction 
of dense gas in the south is significantly lower, which 
should favor HCO"'". That this is not seen implies some 
other excitation mechanism or chemistry is in play in the 
gas of the southwestern arm. From the distribution of 
shock tracers f section |4.4[) . GMC G must be one of the 
clouds most effected by shocks. It is tempting, there- 
fore, to suggest that the southern bar end in Maffei 2 is 
an example of a location of elevated HCN abundances 
in shocks. In the Galaxy, bipolar molecular outflows 
do s how HCN much more strongly elev ated than HCO+ 
(e.g. iBachiller fc Perez GutierrezI [l997l ) . This scenario 
seems less favored when the other dense gas tracers are 
incorporated. Towards GMC G, both the HCO-^/NsH"'- 
and HCO"'"/HNC ratios are anomalously low, while the 
HCN/HNC ratio remains roughly constant relative to the 
rest of the (non-star forming) nuclear disk. This argues 
for a depletion/destruction of HCO"^ rather than an en- 
hancement of HCN. The decrease in HCO"'" intensity here 
may stem from the shock in a more indirect way. The 
passage of a shock alters the chemical state of the molec- 
ular gas, including a dramatic increase in gas-phase wa- 
ter, w hich is effective at destroying HCO"*" (jBergin. et al.l 
|1998[ ). The main difficulty with this mechanism is that 
H2O is also expected to destroy N2H"*", which is not ob- 
served. 

An alternate and perhaps simpler explanation for the 
elevated HCN/HCO"'" ratio and the brightness of high 
density tracers like HNCO and N2H"'" at the southern 
bar end is that while the dense gas fraction decreases 
along the southern bar arm, the density of the dense 
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component does not. In this context the correlation be- 
tween the dense gas tracers and CO/HCN towards the 
star forming ring would imply that the dense gas fraction 
and dense gas density are coupled there. The difference 
between the two environments could be a result of the 
different dynamical environments at each location. To- 
wards the southern bar end one can envision a situation 
where the shocks (section |4AT]) produce a low filling fac- 
tor, high density component that explains the behavior 
of both ratios. The fact that the HNCO/CH3OH ratio 
remains high gives additional direct evidence that den- 
sities of the shocked gas are quite high towards GMC G 
(section HX^ . 

In summary, while the CO/HCN ratio suggests a 
rapidly d ecreasing dense gas fraction along the arms 
(|MTH08f ). The HCN/HCO+ ratio suggests that the ac- 
tual density does not fall significantly. Clearly not all 
the properties of the dense component are traced in the 
CO/HCN ratio. This further emphasizes the fact that 
the CO/HCN measures the fraction of the gas that is 
dense, not necessarily the actual density of the gas. Ra- 
tios that sample only dense gas, such as HCN/HCO'^, 
are needed to determine the physical conditions of the 
dense component. 

4.1.2. HCN/HNC: Proposed Tracer of Warm Gas 

The HCN/HNC line ratio is another dense gas ratio of 
increasing popularity, in which high values indicate hot 
gas and/or PDR conditions. Essentially identical Ay, 
collision coefficients, upper energy states and frequen- 
cies make the ratio particularly insensitive to changing 
excitation. Dissociative recombination of HCNH+ domi- 
nates the formation of both HCN and HNC in cool, dense 
clouds, producing both species in equal amounts (e.g. 
iHirota et al.11199 8'). So in these environments the abun- 
dances of the two species are locked together. Addition- 
ally both transitions likely have significant opacities. It 
is therefore not surprising that this ratio remains near 
unity over a wide range of conditions. When gas tem- 
peratures get large HNC can be converted to HCN via 
HNC -f H ^ HCN -t- H, while in diffuse gas / PDR-like 
environments, the reaction N +CII2 H CN -I- N be- 
comes a major additiona l HCN source (e.g. iSchilke et al.l 
[l99l iTurner et al.llT997l) . 

Fig. [3 (Table HI displays the HCN/HNC intensity ra- 
tio map. Indeed the ratio remains fairly constant across 
the nucleus with values of ^1.1 ±0.3. Variations in the 
ratio are only at the 30 % level over most of the nucleus, 
with the highest values on the star forming eastern ring 
and south of GMC F, which suggests mild enhancement 
by PDRs. Ratios significantly below unity on these small 
spatial scales are not seen anywhere across the nucleus. 
Hence there is no evidence for anomalous enhancement 
of HNC relative to HCN t hat has been found in some 
ULIRGs (|Aalto et al.ll2002[ ). In IC 342, the HCN/HNC 
intensity ratio is also fairly constant , with w eak enhance- 
ment seen towards the PDR GMC (|MT05l ). 

In summary, at least in Maffei 2 and its less active 
cousin, IC 342, the HCN/HNC line ratio has limited di- 
agnostic power. Its diagnostic power becomes stronger 
for molecular gas with extreme temperatures and den- 
sity, such as LIRGs and ULIRGs, where the chemistry is 
driven away from dissociative recombination of HCNH+ . 



4.1.3. HCO+/N2H+: Tracing the Dense Gas 

Both HCN and HCO+ (and to a lesser degree HC3N) 
are observed to maintain large abundances in PDRs, 
shocks and other energetic regions. By contrast, in the 
Galaxy N2H+ is observed to avoid PDR regions, main- 
taining high abundan ces only in dense, qu iescent, dark 
cloud- like conditions (jWomack et al.l ri992j ). In dense, 
quiescent gas both N2H+ and HCO"'" form from proton 
transfer reactions of H;^ with N2 and CO, respectively: 

H+ + N2 ^ N2H+ + H2, (2) 

H+ + CO HCO+ + H2. (3) 

Given the reasonable assumption that the rate 
coefficients for these tw o reactions, } ^\ a nd Ajj], 
are approximately equal ([Woodall et al.l [20071) , then 
X(N2H+)/X(HC0+) ~ X(N2)/X(C0). However, only 
HCO^ maintains a relatively rapid formation route in 
diffuse or PDR gas (C+ -h OH ^ C0+ + H followed by 
H2 + C0+ HCO+ -I- H). No corre sponding route ex- 
ists for N2H+ in PDRs (Turner"1995b). This explains the 
observed behavior of the HCO+/N2H+ fine ratio (Fig. [7] 
and Table [1]) . Towards the PDR regions associated with 
the active star formation sites, HCO^ can maintain ap- 
preciable abundances while N2H+ is readily destroyed. 
The fact that this ratio, which reflects dense gas prop- 
erties, is modified by the star burst PDRs implies that 
towards this region (GMC E and the eastern ring) PDRs 
penetrate the dense gas. The strong correlation between 
high HCO^ /N2H^ ratios and PDR tracers like C2H, sug- 
gests that this ratio may he a powerful new probe 0/ dense 
PDRs in external galaxies. 

Along the arms, X(N2H+)/Ar(HC0+) abundances are 
typically ^0.1-0.15, with values reaching up to ^0.33 to- 
wards the southern bar end. If a standard CO abun- 
dance of X{GO) ~ 10~^ is adopted, Ar(N2) abundances 
of 1 — 3 X 10~^ are implie d. If metallicities are solar (e.g. 
lAnders fc Grevess3ll989( ) and all N is in the form of N2 , 
then X(N2) ~ 4 x lO''^. Evidently, in the more qui- 
escent dense clouds of the nucleus of Maffei 2, at least 
a quarter of all nitrogen is in molecular form. Using 
the above X(N2) estimation together with Table [1] gives 
X(N2)/A:(N2H+) ~ 2 X 10^ for these regions. 

Following the analysis ofiMTOS, X(N2)/A:(N2H+) pro- 
vides indirect constraints on the cosmic ray ionization 
rate, C- Adopting A:(N2)/X(N2H+) - 2 x 10^ and 
the extremely conservative lower limit, X{e^) = X{ion) 
> A:(N 2H+) + X(HCO+) > 6 - 9 X 10-^ Figure 9 of 
IMT05I indicates C/nH2 > 10~^^ s^^ cm^. Less conserva- 
tively, for X{e~) ~ 10~^ typical of Galactic star forming 
regions fe.g. iCaselh eFal] [T998[) . C/uh^ > lO'^"-^ s'^ 
cm"^. Therefore, ^ is at least the Galactic value if the 
density of the gas that the N2H''' and HCO^ emission 
originates in is > 10"* cm^"^. These constraints are 
comparable to though somewhat stronger than observed 
for IC 342. 

4.1.4. The GS(2-l)/HCN(l-0) and CS(2-l)/(f^ S(2-l) 

Ratios 

iKuno et al.l (|2008[ ) has pubhshed a map of CS(2-1) to- 
wards the nucleus of Maffei 2 with similar spatial reso- 
lution to the data presented here. We make use of this 
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map to compute the CS/C^''S and CS/HCN ratios to- 
wards the GMCs. The CS/HCN intensity ratio is fairly 
constant across the nucleus except towards the southern 
arm (GMC G), where it is a factor of two lower. If strong 
chemical evolution has not altered the '^^S/'^'^S ratio, the 
C^'^S emission implies the regions of high CS opacity are 
along the southern arm and not at the CS peak near 
the massive star formation. This is unexpected. For a 
^^S/^^'S ratio of 23, LTE hue ratios imply CS opacities 
of '-~^6 along the southern arm. Comparing the observed 
peak temperature of C'^^S gas to optically thin emission 
with Tea: = 10 K yields a filling factor, /(C^^S) ~ 0.018. 
Given the differences in beam size thi s filling f actor is in- 
termediate between ^^CO and HCN (IMTH08D . Towards 
the GMCs D+E the larger CS/C^'^S ratio implies that 
CS opacities here must be less than 1.4. Since the inten- 
sity of the main isotopologue peaks towards this region, 
gas temperatures or filling factors are much larger here. 
So in order to explain the morphology of both CS and 
C^'^S, it is required that CS be warmer and more dif- 
fuse towards GMCs D-l-E when compared to the GMC 
G. In summary CS appears compatible with the inter- 
pretation based on HCN/HCO"'" that the southern bar 
end has high density but low dense gas fraction. 

4.1.5. HC-iN(10-9) 

IIC3N has a critical density that is comparable or 
slightly higher than HCN but very likely has low opacity 
due to its lower abundance. Moreover its levels are more 
closely spaced and so population is distributed over a 
wider range of levels, making any one level more sensitive 
to changes in excitation. HC3N originates primarily from 
the arm-ring intersection regions similar to HCN and 
HCO+, although perhaps peaks slightly farther from the 
center of the galaxy than does HCN. Using the kinemat- 
ics to further resolve the emission suggests that HC3N 
is tracing the inner portion of the arm rather than the 
Xi — xi intersection region. Evidently, densities increase 
along the arms as the central ring is approached. Unlike 
HCN, HC3N weakens just before the strong star forma- 
tion is reached at GMC D-l-E and just north of GMC 
F. Potentially this results from warmer temperatures de- 
populating the J = 10 level towards the starburst. 

4.2. PDR Chemistry and the Nuclear Starburst 

Ionizing radiation from massive stars impinging on 
nearby molecular clouds changes their chemical compo- 
sition. Large abundances of ionized carbon in PDRs 
drive a rich hydrocarbon chemistry on cloud surfaces. 
One of the simplest species from such regions is C2H. 
C2H is formed from rapid dissociative recombination 
of C2Hj" c oming from reactions be tween C+ + CH3 
(|Wootten et a l. 1980; S ternberg fc Dal garno 1995). C2H 
chemistry is particularly closely tied to ionized carbon 
and is an excellent probe of these PDRs. Its usefu lness as 
an extragalactic PDR tracer was demonstrated i nlMTOSl 

After accounting for opacity effects (see H3.3\ for 
C2H(l-0; h — ^) for an optically thin C2H transition) 
C2H pervades much of the molecular clouds in the north- 
ern star forming region. Lower level C2H emission ex- 
tends over the rest of the nucleus but avoids the cloud 
cores, instead favoring the side of the dense cloud facing 
star formation (traced here by 2 cm radio continuum) 



(see Fig. [S]). This is particularly clearly seen west of 
GMC F and north east of GMC G. 

As discussed in i i4.1.31 the dramatic increase in the 
HCO+/N2H+ line ratio across the central ring shows 
not only that PDRs are prevalent near the young, star 
forming region but that the PDRs penetrate much of 
the dense gas there. Given the distinct region of al- 
tered chemistry, constraints can be placed on the range 
of influence of the current burst of star formation. For 
the northern star forming region, the region of elevated 
HC0'^/N2H+ ratio is 80 pc diameter in extent. We as- 
sume that the volume of the region emitting strong PDR 
emission reflects the volume of molecular gas engulfed 
within HII regions that have expanded into the ambient 
medium due to over-pressure. This should be reasonable 
on GMC scales. For an ionized gas sound speed of 10 
km s^-'^, N^yc = 5 X 10^^ s~^ and a uniform density of 
^ 10^-'^^ cm"3 typical of the central ring (|MTH08l ). 
a standard expanding HII region would reach a radius 
of 40 pc in ~7 Myr. While this analysis is sensitive to 
the exact density distribution, it is interesting that this 
timescale for dynamical evolution of GMCs due to pres- 
sure is an order of magnitude longer than characteristic 
chemical timescales and consistent with the timescale of 
massive star evolution. 



4.3. The Molecular Outflow along the Minor Axis in 
PDR Gas 

C2H emission appears to trace a molecular outflow 
along the minor axis (Fig. The northwest exten- 

sion of this bipolar flow is redshifted by ^ 15 km s~^ 
and the southeast blob blueshifted by ~ —15 km 
with respect to the velocity of the corresponding in-plane 
material. These velocity shifts have the correct sign for 
outflow given Maffei 2's geometry (northern arm being 
nearer to us, and the arms trailing). There is tenta- 
tive evide nce for th e flow in HNC and SiO (and possibly 
CO(l-0): [MTH08l ) but the molecular outflow is not ob- 
vious in any of the other species. This is understandable 
given the strength of C2H in the outflow reflects its rela- 
tively low critical density and its ability to maintain high 
abundances in the presence of ionizing radiation. 

The molecular outflow in Maffei 2 does not appear to 
follow the distribution of ionized gas traced by 2 cm. 
This is different from wha t is seen in the muc h more pro- 
nounced outflow in M82 (jWalter et al.ll2002t ). The 2 cm 
radio continuum emission appears to trace t he walls of 
Maffei 2's outflow (Fig. l9l I Turner fc Hoi HOOl . The C2H 
emission must trace entrained molecular material being 
driven out directly. The explanation for the 2 cm emis- 
sion along the edges remains unclear but may be due to 
a path length effect of a cone of ionized gas surrounding 
the molecular flow. Since evidence is seen for a small 
minor axis outflow it is possible to obtain a second esti- 
mate of the characteristic age by assuming that the gas 
traced in C2H is free to flow ballistically along the mi- 
nor axis. The extent of the outflow is ±200 pc. For an 
inclination angle of 67° and the measured line of sight 
velocity, the time of flight corresponds to 6 Myr. This 
is a very young outflow. The agreement with the char- 
acteristic in-plane timescale is excellent, suggesting that 
the current starburst may have an age of r ~ 6 - 7 Myr. 
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4.4. Shocks and Gas- Grain Chemistry: GH^OH, HNCO 

and SiO 

Elevated gas-phase abundances of species whose forma- 
tion is sensitive to the grain processing, such as CH3OH, 
HNCO, or SiO, are impor tant p robes of hot core gas 
and shocks. As noted in IMT05I . and illustrated even 
more clearly here, these species do not trace hot cores. 
Typical Galactic hot cores, even significant collections 
of them cannot produce the required emission necessary 
to explain the intensities on these large scales. More- 
over, we detect no spatial connection between CH3OH 
and HNCO, and star formation, if anything, we see an 
anticorrelation. We can, therefore, safely conclude that 
the above three species are not tracing hot core emis- 
sion or outflows from young stellar objects but must 
originate in large scale shocks. Inter-comparis qns can 
be used to identify the resp o nsible shock source (IMTOa 
Requena- Torres et al.' 1 20061 : iRodriguez-Fernandez et all 
2006; Martin et al. 2008|) and place constraints on shock 
properties such as shock strength, type of shock and lib- 
erated grain / mantle composition (.Usero et al. ,2006). 

In this dataset three shock tracers, CH3OH, HNCO 
and SiO, are mapped. Each is sensitive to somewhat 
different shock properties. It is expected that elevated 
gas-phase SiO results from release of Si from destruction 
of silicate grain cores (e.g. Martfn-Pintado et al. 1992). 
That the grain core must be disrupted requires 'strong' 
shocks capable of sputtering grain cores. High observed 
abundances of gas- phase CH3OH likely come from direct 
liberation of CH3OH that has formed on grain surfaces 
by complete hydroge nation of CO fe.g. lMillar et al.iri991l : 
IBachiller et al.|[T995() . Liberation of mantles does not re- 
quire shocks as strong as those that effect the grain core. 
Thus the relative enhancement of SiO versus CH3OH 
will depend on the strength of the shock. Grain core de- 
struction requires shock velocities of s^25 km s~^ whereas 
shock velocities of only > 10-15 km s~^ are required 
for s ignificantly elevated gas-phase CH3 OH abundances 
fe.g. ICaselh et al.llT997l : IBergin. et al.lll998h . It is rea- 
sonably well establis hed that HNCO is a mantle species 
along with CH3OH (|MT05t iRodn'guez-Fernandez et al.l 
l2010f) . But its formation mechanism remains unclear. 
The most commonly proposed mechanism is partial hy- 
drogenation of OCN on grain surfaces followed by direct 
liberation (e.g. Hs&Leg_awa_&J3£xbst, JJiQS,). Other sug- 
gested mechanisms include acid-b ase reactions between 
CO and NH3 and UV ph otolysis (jSchutte k, Greenbergj 
[1997I: IKeane et al.l [200I . Thus if HNCO forms on 
grain mantles and is liberated directly then its unsat- 
urated nature implies that it originates in grain man- 
tles with diff erent chemical c ompos ition from CH3OH 
(however see iTideswell et al.l (|2010D for a possible for- 
mation mechanism in more saturated ices). Unsat- 
urated ices are favored when gas accretes in atomic 
hydrogen-poor environment whereas sat urated ices fa- 
vor a tomic hydrogen-rich accretion (e.g. Cha rnlev et al.l 
[19971 [Wata nabc fc Kouchi.2002. ). Hydro gen-poor condi- 
tions reflect accretion at high densities when almost all 
hydrogen is in molecular form. If it is assumed that the 
shock simply injects molecules into the gas phase, the 
HNCO/CH3OH abundance ratio is then sensitive to gas 
physical conditions during ice mantle formation. 

Chemistry is not the only effect that influences the 



relative abundance of these three species. Excitation is 
also important. The excitation requirements of SiO(2- 
1), CH30H(2fc - Ifc) and HNCO(5o5 - ^m) transitions 
differ (see Table [J). The critical density of the HNCO 
transition is an order of magnitude higher than SiO and 
CH3OH. Furthermore the energy of the upper state in 
HNCO(5o5 — 4o4) is more than twice that of SiO and 
higher than two of the three main transitions that make 
up the CH3OH line. Therefore, CH3OH and SiO lines 
will be favored over HNCO in lower density gas. Also 
in the presence of ionizing radiation each species disso- 
ciates at different rates. HNCO is photodissociatcd at 
a rate twice that of CH3OH and ■-^SO times that of SiO 
(jRoberge et al.lll991l: [Sternberg fc Dalgarndll995[ ). 

Thus considering the distribution of CH3OH, HNCO, 
and SiO can give insights into both the chemical and 
physical conditions present in the shocked gas and po- 
tentially tell us about nature of the shocks themselves. 
Below we discuss the overall structure of the shock trac- 
ers in the context of nuclear bar models, then we analyze 
the subtle differences between each shock tracer to fur- 
ther constrain shock properties. 

4.4.1. Overall Shock Tracer Morphology and Bar Structure 

In Maffei 2 the outer bar arms dominate the overall 
distribution of the shock tracers CH3OH, HNCO, and 
SiO, whereas emission from the star forming ring is min- 
imal. This is very different from the morphology of CO, 
HCN and the PDR tracers. The bar ends are especially 
bright in the shock tracers. (The faintness of shock trac- 
ers along the northern arm is presumably due to the fact 
that there is much less molecular gas there, particularly 
at what would be the location of the northern bar end.) 
It is perhaps surprising that it is the bar ends and not 
the inner arms that have the largest fraction of the gas 
influenced by shocks. This provides constraints on the 
nature of gas flows in Maffei 2. 

Different models exist for describing the structure of 
gas bars, each predicting different shock strengths and 
locations. The models fall into two broad classes, orbit- 
crowding models and hydrodynamical flows. Orbit- 
crowding models treat the gas as distinct clouds orbiting 
on essentially ballistic orbits set up by th e potential (e.g. 
ICombes fc GCTh][l985l : [Bvrd et al.l[T99l . The arm mor- 
phology of the bar is then determined by where orbits 
converge. In hydrodynamical models the gas is treated 
as a fluid that responds to the stellar potential and is 
susceptible to hydrodynamical forces such as pressure 
and viscosity. In such models, gas flows along stream- 
lines with shocks occurring where the streamli nes col- 
lide (e.g. [Athanaiiouii [Tl9l [EiiEiEin [l^^ Both 
models predict similar morphologic al distributions bu t 
rather distinct shock properties (e.g. iRegan et al.l ll999\ 
In orbit-crowding models the bar arms result from molec- 
ular clouds moving along converging orbits. In such mod- 
els only weak shocks are expected because the bulk of the 
velocities are parallel to both clouds (with the possible 
exception of arm / nuclear ring intersections) . Moreover 
such models do not isolate the bar ends as locations of 
strongly convergent orbits and so it is not expected that 
these location would be sites of vigorous shock activity. 

Hydrodynamical models, on the other hand, show 
gas flow streamlines that change directions dramatically. 
The places where flows converge are accompanied by 
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strong changes in velocity, and hence strong shocks. In 
hydro models, gas drifts slowly inward along bar-aligned 
xi stellar orbits (or in the case of Maffei 2, from radial 
inflow along the outer bar) until the gas reaches the in- 
nermost, cuspy xi orbit (of the nuclear bar). The cusps 
of these orbits lie at the bar ends. There the gas is com- 
pressed, shocks and plunges nearly radially inward along 
the bar arm until it reaches the much further in, perpen- 
dicular to the bar X2 orbitals, where it is again shocked, 
loses angular momentum and transitions onto the X2 or- 
bitals. Some high angular momentum gas 'sprays' off 
the inner X2 orbit s, swinging out and colliding with the 
arm gas fl ow (e g. IAthanassoulalll992t iPiner et al.lll995l: 
iKim et aLll2011[ ). The amount, location and relative ve- 
locity of gas in the spray region depends on bar param- 
eters. 

In summary, the strength of the emission from the 
shock tracers CH3OH, HNCO, and SiO along the bright 
bar ends favors the hydrodynamical bar models, as orbit- 
crowding models do not predict strong cloud collisions at 
the bar end whereas they are sites of shocks in hydrody- 
namical models. 

4.4.2. Relative Shock Tracer Morphologies and Bar 
Properties 

What about changes in shock properties across the 
nuclear bar? Fig. [TO] displays the HNCO/CH3OH, 
SiO/CHgOH and SiO/HNCO peak temperature ratios 
smoothed to 8 and sampled at the GMCs. Given the 
significant difference in the native resolution of the three 
datasets and the fact that they were obtained with two 
different interferometers, only the clearest trends are dis- 
cussed. 

As demonstrated by the principal component analysis, 
the morphologies of the three species are strongly cor- 
related, implying that we have a robust measure of the 
locations of shocked gas. So similar are the maps that 
only very subtle differences between HNCO, CH3OH and 
SiO appear, suggesting that shock conditions are uniform 
across the nucleus. With the possible exception of the 
uncertain SiO ratios towards GMC B, the arm clouds (C, 
D, F, G) have effectively constant HNCO/CH3OH (- 
1.1 ± 0.2), SiO/CHgOH (- 0.20 ± 0.02) and SiO/HNCO 
(- 0.19 ± 0.03) ratios. Only towards GMC E and the off 
arm spray regions (GMC H and just south of the west- 
ern ring) are different ratios observed. Towards GMC 
E (the northern star forming region), HNCO/CH3OH is 
lower and the SiO ratios are higher relative to the rest 
of the arm. In the off arm regions HNCO/CII3OH ratios 
are depressed while the SiO/CIl30H ratio is elevated and 
SiO/HNCO are elevated (spray region) or flat (GMC H). 

From these ratios, three main points are concluded. 1) 
While GMC G (the bar end) is the brightest source in 
the nuclear bar, indicating that the largest fraction of 
gas along this line of sight is shocked, evidently shock 
strengths are not enhanced here. SiO to mantle species 
have the same basic ratio at the bar end as it has over 
the rest of the bar. 

2) Towards GMC E, the decreased intensity of HNCO 
relative to CH3OH and SiO is consistent with what is 
expected from photodissociation. From section 14.21 we 
know that the dense gas here is permeated by UV pho- 
tons for the young star formation. The shock tracer ra- 
tios further conflrm this. Because of this, the absence 



(or faintness) of shock tracers towards the central ring 
and in particular the northern xi — X2 orbit intersection 
cannot be taken as evidence of absence of shocks here. 

3) The spray region upstream (counterclockwise) from 
the main bar arms favor CH3OH and to a lesser degree 
SiO relative to HNCO. It is possible that because gas 
is more diffuse in this region grain mantles are H-richer 
than in the denser arm. However, we favor a simpler 
explanation associated with changing excitation. The 
spray regions have lower density hence the higher critical 
density HNCO transition is expected to be weaker rela- 
tive to the lower critical density, CH3OH and SiO tran- 
sitions. That the Sqs — 4o4 transition of HNCO preferen- 
tially favors the dense clouds relative to more extended 
CH3OH explains why larger beam single-dish multi-line 
observations derive an excitation temperature for HNCO 
that is nearly a factor of ^ 2 larger than derived for the 
CH3OH transitions in Maffei 2 (Hiittcme ister et al.lll997l : 
iMartfn et al.l [20091) . Note that the excitation tempera- 
ture of both species is quite low (Te^ <15 K), implying 
that the shock regions remain cool or the densities are 
low enough for these transitions to be sub-thermal. 

In summary, we conclude that there is little evidence 
for strongly varying grain chemistry across the shocked 
regions on 8 scales. The small level of variation seen 
is consistent with changing physical conditions (radiation 
field and gas density). 

4.5. The Chemistry of Spiral Nuclei: Maffei 2 vs. IC 

342 

The overall nuclear structures of Maffei 2 and IC 342 
are similar. Both have extended bar-like arms that termi- 
nate on a central ring that is actively forming stars. The 
molecular gas morphol ogy of IC 342 is con sistent with 
that of a gaseous bar ([Turner fc HurtI [19921 ) . However, 
because of the very small central ring and ambiguous 
kinematics, other conclusions including unbarred spiral 
arm s coupled with massive star feedback remain feasi- 
ble ([Schinnerer et al.l 120081 ). On the other hand Maffei 
2 is clearly a nuclear bar bas e d on both i ts morphology 
and kinematics fe.g. lMTH0"8l : iKuno et"ar.l,2008 '). Given 
the potential dynamical similarities of the two nuclei a 
comparison of the chemistry of the two can provide im- 
portant insights on the genera hzabili ty of the the overall 
chemical picture developed in lMTOSl 

The statistically robust principal component analysis 
demonstrates that the chemistry of the two nuclei are 
remarkably consistent. Maffei 2 can be deflned chemi- 
cally as a highly inclined copy of IC 342 (Fig. [TTl . In 
both systems the first principal component characterizes 
the density-weighted column density, and peaks along 
the inner arms and central ring. The second principal 
component separates the main chemical forces responsi- 
ble for controlling observed abundances, namely shocks 
versus PDRs. PDR species dominate close to the young 
star forming region while the shock tracers peak along 
the arm, with preference for the cuspy bar ends. As a 
result, the bar driv en star formation / chemical scenario 
outlined in IMT05I must apply also to Maffei 2's nuclear 
bar. 

However the chemical picture is not precisely identical, 
and these differences give added insight into the relative 
differences in evolutionary phase of the two nuclei. The 
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main difference between the two galaxies is in the distri- 
bution of N2H+. In IC 342 the PDR tracers are confined, 
appearing to be primarily cloud surface features, lining 
the inner ring. Significant N2H+ (and HNCO) seen near 
the young star forming region in IC 342 implies that its 
main star forming region influences the dense gas only 
locally, on scales of 50 pc or less. Whereas in Maffei 2, 
the N2H^ /HCO^ ratio and photodissociation of HNCO 
clearly demonstrate that the star formation has begun to 
dramatically alter the chemical state of the dense molec- 
ular clouds in its immediate environment, not just the 
surfaces. 

Evidently the central molecular ring in Maffei 2 is 
structurally different (lower clump density) and that to- 
wards CMC E star formation is somewhat more evolved 
than the main burst in IC 342. This is consi stent with the 
steeper radio spectrum towards CMC E (jTurner fc Hoi 
[T994i: .Tsai et al....200Ql . The crude estimates of timescale 
required to alter the chemistry of the dense component 
is ~ 7 Myrs. This is a reasonable timescale compared to 
stellar evolution of the burst in Maffei 2 and the kinemat- 
ics of the outflow. However, this timescale is still very 
short; short enough that even this 'evolved' part of the 
burst is in its early stages. The star formation extending 
to the north and south of CMC E is likely as young as 
the embedded star formation in IC 342. 

In IC 342 the peak temperature ratios between SiO, 
CH3OH and HNCO, are similar to those seen in Maf- 
fei 2. The PDR CMC in IC 342 has slightly lower 
HNCO/CH3OH and elevated SiO/HNCO ratios. This 
is exactly the behavior exhibited by the FDR-modified 
CMC E in Maffei 2, lending strong support for the in- 
fluence of photodissociation on shock tracer brightness. 
The shock chemistry dominated cloud in IC 342 has 
slightly elevated HNCO relative to the other species, 
making its set of ratios a match within the uncertain- 
ties of CMC G. This agreement suggests that both the 
grain composition and shock properties are the same 
for these structurally similar locales. However, since 
brightnesses are sensitive to radiative transfer effects, 
one must be careful over-interpreting ratios if gas excita- 
tion is changing rapidly. Maffei 2 and IC 342 do have 
rather similar CO excitat i on cond itions (| Meier et al.l 
[2OOOI: IMeier fc Turned[20Qll: IMTH08D . 

In summary, we find that when the dynamical environ- 
ments are similar, as they are in IC 342 and Maffei 2, a 
consistent chemical picture begins to emerge. Gas chem- 
istry is confirmed as a powerful method for isolating the 
physical mechanisms that influence gas conditions in spi- 
ral nuclei. In future papers in this series, we will broaden 
the range of galactic environments studied to assess how 
far chemistry can be pushed as a diagnostic of sub-beam 
physics. 

5. CONCLUSIONS 

This paper presents 2" - 10" imaging of transitions 
from nine separate molecules across the nuclear bar in 
Maffei 2 obtained with the BIMA and OVRO millimeter 
interferometers. The results are compared with the first 
paper in the series on IC 342 to test the generalizability 
of the chemistry discovered there and to constrain dif- 
ferences in structure between the two nuclei. Our main 
conclusions are as follows: 



1. Morphologies of the observed transitions exhibit 
both dramatic (e.g., C2H vs. CH3OH) and subtle 
(e.g., HCN vs. HCO^) spatial differences indica- 
tive of the varying dominance of shocks vs. radia- 
tion fields as drivers of the chemistry in individual 
clouds. 

2. A principal component analysis reveals that HCN, 
HNC, HCO+ and 3 mm continuum are all ex- 
tremely tightly correlated, indicating a direct con- 
nection to massive star formation. 

3. CH3OH and HNCO are tightly correlated with 
each other and are strongly anti-correlated with 
the tracers of massive star formation; based on 
their spatial distribution, it appears that these two 
molecules are formed by liberation of grain mantles 
in the presence of bar shocks. 

4. Variations in the ratios amongst the popular A ^ 3 
mm dense gas tracers, HCN, HCO+, HNC, CS and 
HC3N, are at or less than a factor of three, unlike 
the order of magnitude variations exhibited by the 
chemistry-sensitive molecules. This indicates that 
variations in excitation in the cool, dense gas are 
much less dramatic than variations in dense gas 
chemical conditions. 

5. Localized variations in the ratios HCN/HCO+, 
HCN/HNC, HCN/CS and HCN/CO can be ex- 
plained by variations in density. The dense gas 
ratios indicate that the highest density gas is found 
at the xi —X2 orbit intersections, the lowest density 
gas is along the non star-forming western portion 
of the ring, and intermediate density gas is found 
near the star formation. Although high density gas 
is found throughout the full extent of the nuclear 
bar arms, the dense gas fraction falls at larger radii 
from the nucleus. 

6. There is no evidence for X-ray dissociation regions 
or AGN infiuences on the gas in the center of Maf- 
fei 2, specifically, as traced by the HCN/HCO+ ra- 
tio. 

7. The strong correlation between high HCO+/N2H+ 
ratios and the PDR tracer molecule, C2H, suggests 
that this ratio may be a powerful new probe of 
dense PDRs in external galaxies. Given the high 
critical density of both HCO+ and N2H+ m com- 
parison to C2H, the elevated ratio demonstrates 
that PDR gas towards the strong star forming re- 
gions permeate the dense clumps. When combined 
with the morphology of C2H this ratio may provide 
useful insights into the relative density of PDR gas. 

8. C2H emission reveals a molecular outflow along 
the minor axis that has entrained molecular gas. 
Wrapped around the outside of this outflow is ion- 
ized gas traced in 2 cm radio continuum. Assuming 
ballistic motion, the morphology and kinematics 
are consistent with a very young age of 6-7 Myrs 
for the outflow. 

9. The bar induced orbital shocks in Maffei 2 are 
traced in three species, HNCO, CH3OH and SiO. 
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All three species support the conclusion that the 
nuclear bar end is the site of the largest fraction 
of shocked gas. That this location exhibits such 
strong shock properties argues in favor of hydrody- 
namical bar models. 
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Fig. 1. — The basic structure of MafEei 2's nucleus. Left) The high resolution ^^CO(l-O) integrated intensity map from lMTHOl with the 
location of the bar arms shown (in blue: online edition). Contours are 5.0 Jy bm~^ km s~^ X 1, 2, 4, 8, 1 6. A green cross marks the 
dynamical center. Right) The same CO contours overlaid on a greyscale of the 89 GHz continuum IIMTHOSI) . The greyscale ranges from 
0.5 - 5.0 mjy bm'^. The red crosses represent the location of the labeled GMCs. The location of the central ring and active star formation 
complexes are also noted. 



TABLE 1 

Molecule and Column Density Parameters 



Molecular 


M 


A/B/C 


SuWKgi 


Eu 




^mol / -^mol 


iV(50) b 
]V(10) 


Transition 


(Dby) 


(GHz) 




(K) 


log( cm~^) 


/ cm^ ) 
' Kkms ~^ ^ 




SiO(2 -l;v=l) 


3.098 


■■■ /21.712/ ■•• 


2 


1787 


... / ... 






SiO(2 -1; v = 0) 


3.098 


•■■ /21.712/ ■■• 


2 


6.25 


5.27/5.66 


1.86(12) 


2.94 


C2H{l-0; |-|) 
C2H{l-0; i-l) 


0.8 


■ ■ ■ /43.675/ ■ • • 


1 


4.19 


5.13/6.24 


2.32(13) 


3.01 


0.8 


■ ■ ■ /43.675/ ■ • • 


1 


4.19 


5.13/6.24 


4.64(13) 


3.01 


HCO+(l - 0) 


3.888 


. . . /44.594/ ■ • • 


1 


4.28 


5.15/5.80 


9.51(11) 


3.36 


HNC(1 - 0) 


3.05 


■ ■ ■ /45.332/ ■ • • 


1 


4.35 


6.39/5.87 


1.51(12) 


2.94 


HC3N{10 - 9) 


3.72 


•■■ /4.5491/ ■•• 


10 


24.02 


5.86/6.87 


6.73(12) 


0.723 


N2H+(1 - 0) 


3.40 


• ■ ■ /46.587/ ■ • • 


1 


4.47 


5.63/5.85 


1.16(12) 


3.00 


C34S{2 - 1) 


1.96 


• ■ ■ /24.104/ ■ ■ ■ 


2 


6.94 


5.65/6.17 


4.06(12) 


2.56 


CH30H(2fc - Ifc) 


0.89 


127.484/24.680/23.770 


4,3,4== 


6.98 


4.83/6.39 


1.58(13) 


5.96 


HNCO(5o5 - 4o4) 


1.60 


912.711/11.071/10.911 


5 


15.83 


6.29/5.99 


9.27(12) 


3.15 



Note. — Data from Splatalogue, the JPL Molecular Spectroscopy Catalog ^ickctt ct al. 1998) and references within. Collisional coefficients are 
C^'^S, HC3N HGreen fc ChapmanllToTa) . N2H+ l IGreenlfToTa) . CH3OH l lPottage ct al...20p 4 ,). HNCO jGreenlfToga) . HCO+ jFlowei|[T999h and SiO 
IDavou fc Balancall2006f) . ForllNC and C2H we adopt the collision coefficients of HON llCreen fc ThaddeuJl97l) . 

^ The critical density of the transitions neglecting opacity effects {^^ricr — TT^^ Sind '^n^r based on the formalism of IDickinson et al.l ]197711 ■ with 
[c^] ~ 1.0 X 10"^. 

^ The ratio by which the derived column densities change when the Tex is changed from 10 K — ^ 50 K. 
SuISkQi for the three blended 2i — liE, 2o — IqE, and 2o — I0A+ transitions (eg.. lTurnei1ll99 J) . 



Resolved Chemistry in Maffei 2 



15 




Fig. 2. — Spectra of t he molec ular transitions towards Maffei 2's GMCs. Spectra are summed over an 8 aperture centered on the GMC 
positions indentified bv lMTffOSI . The vertical axis is the flux in the aperture and the x axis is LSR velocity. 
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Fig. 3. — Integrated intensity maps of four molecular transitions in Maffei 2 observed with OVRO. a) Tlie ^'^CO(l-O) in steps of 0.75 Jy 
bm-l kms-l. I IJVLTH08I) b) The HCN(l-O) I IJVLTH08I) insteps of 0.50 Jy bm"! kms"! with a resolution of 2.5" x2.2". c) The HCO+(1-0) 
in steps of 0.50 Jy bm~^ km s~^ for a resolution given in Tablc[T] d) The HNCO(5o5-4o4) in steps of 0.40 Jy bm^^ km s~^. In all planes 
the greyscale is the ^"^00(1-0) integrated intensity. 
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Fig. 1. Integrated intensity maps of the molecular transitions in Mallei 2 observed with BIMA. a) The HNC(l-O) in steps of 1.0 Jy 
bm-i km 8"^ b) The C2H(l-0;3/2-l/2) in steps of 1.1 Jy bm"! km 8"^ c) The C2H(l-0;l/2-l/2) in steps of 1.1 Jy bm"! km 8"^ d) 
The C^*S(2-1) in steps of 1.2 Jy bm~^ km . In all planes the greyscale is the ^^CO(l-O) integrated intensity. 
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Fig. 4. — Integrated intensity maps of the molecular transitions in MafTei 2 observed with BIMA, continued, e) The HC3N(10— 9) 
integrated intensity in steps of levels 0.60 Jy bm"! km s'^ f) The N2H+(l-0) in steps of 0.75 Jy bm-^ km s'^ g) The SiO(2-l; v=0) 
in steps of 0.40 Jy bm~i km s~^. h) The CH30H(2j,-lfc) in steps of 1.5 Jy bm"'^ km . In all planes the greyscale is the ^*CO(1-0) 
integrated intensity. 
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Fig. 5. — Maps of the first three principal components of the molecular distribution. 



TABLE 2 
Observational Data 



Transition 
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Beam 


K/Jy 


Noise 




(MMYY) 




(GHz) 
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{km s~^) 


( X 




(mJy bm~^) 



BIMA; 


















SiO(2-l;v=0)'» 
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B,C,D 
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10 
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11 
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8.1 X 6.5;-86 
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15 
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9.5 


HCO+(1-0)'' 


0199-0399 


L,H,UH 
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13.44 


2.5 X 2.1;-24 


30.5 


6.7 


HNCO(5oB-4o4)'= 


1098-0199 


L,H 


109.906 


240-430 


10.92 
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= 02''41'"55^0 S 
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TABLE 3 

Measured Intensities, Linewidths & Centroids 



GMC 


Transition 




Ti 


Vo 


Av 






(K km s-i) 


{mK) 


{km s~^) 


{km s~^) 



A 

(02:41:57.17) 
(59:36:42.4) 



(02:41:55.75) 
(59:36:29.6) 



(02:41:55.48) 
(59:36:24.8) 



D 

(02:41:55.14) 
(59:36:20.7) 



(02:41:55.08) 
(59:36:18.4) 



F 

(02:41:54.86) 
(59:36:09.6) 



SiO(l-0:v=l) 
SiO(l-0;v=0) 
C2H(|-i) 
C2H(I-I) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
C34S(2-1) 
CH30H(2fc-lfc) 
HNCO(5o6-4o4) 
SiO(l-0;v=l) 
SiO(l-0;v=0) 
C2H(|-i) 
C2H(|-I) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fc-lfc) 
C34S(2-1) 

HNCO(5o5-4o4) 
SiO(l-0;v=l) 
SiO(l-0;v=0) 
C2H(|-i) 
C2H(|-|) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fc-lfc) 
C34S(2-1) 

HNCO(5o5-4o4) 
SiO(l-0;v=l) 
SiO(l-0;v=0) 

C2H(|-I) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fe-lfe) 

e^4s(2-i) 

HNCO(5()5-4()4) 
SiO(l-0;v=l) 
SiO(l-0;v=0) 
C2H(|-i) 
C2H(I-I) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fe-lfe) 

0343(2-1) 
HNCO(5o5-4o4) 
SiO(l-0;v=l) 
SiO(l-0;v=0) 
C2H(|-i) 
C2H(|-|) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fc-lfe) 

C-^4S(2-1) 
HNCO(5o5-4o4) 



<19 
<1.0 
<3.3 
<3.3 
16±7 
<2.9 
<1.0 
<1.9 
<2.2 
<3.2 

<19 
<1.0 
<3.3 
<3.3 
61±7 

9.1±1 
<1.0 

5.1±1 

10±2 

<2.2 

^4.8 

<19 

<1.0 

11±2 

3.9±2 
82±7 
24±1 
2.8±0.5 

9.6±1 

5.0±2 
<2.2 

;S4.8 

<19 
2.5±0.5 

17±2 

11±2 
180±7 
41±1 

3.2±0.5 
9.4±1 
5.5±2 
<2.2 
<4.8 
<19 
~0.77 
19±2 
13±2 
130±7 
44±1 

2.1±0.5 
6.6±1 
6.7±2 
<2.2 
<4.8 
<19 
<1.0 
9.4±2 
<3.3 
85±7 
36±1 

4.1±0.5 
15±1 
15±2 
5. Oil 
22±2 



<88 
<23 
<34 
<34 
270±38 
45±16 
<26 
<30 
<28 
46±19 

<88 
35±12 

<34 

<34 
470±38 
150±16 

<26 
60±15 
53±19 

<28 
50±29 

<88 
24±12 
100±17 
48±17 
520±38 
200±16 
42±13 
110±15 
100±19 

<28 
130±29 

<88 

<23 
110±17 
80±17 
580±38 
250±16 
50±13 
110±15 
95±19 

<28 
86±29 

<88 
22±12 
97±17 
77±17 
520±38 
220±16 
39±13 
85±15 
70±19 

<28 

<58 

<88 
35±12 
63±17 

~46 
450±38 
250±16 
31±13 
150±15 
190±19 

~30 
190±29 



-34±5 
-lOOitll 



-100±6 



-67±13 



-84±3 
-134±3 

-98±6 
-87±8 

-112±20 

-80±20 
-88±6 
-76±17 
-78±3 
-120±4 
-88±7 
-88±5 
-99±4 

-89±5 



-82±9 
-77±17 
-74±3 

-109±4 
-80±6 
-87±4 

-101±5 

-90±11 

-2.0±26 
-73±12 
-72±16 
-73±4 
-104±4 
-80±8 
-83±5 
-94±8 



17±8 
10±12 

-6.9±4 
-0.9±3 
0.1±13 
5.4±6 
-43±3 
-46±26 
-1.6±3 



21±10 
120±25 



45±15 



55±30 



68±7 
65±7 

50±10 
48±20 

110±50 

83±47 
93±15 
90±40 

83±6 
120±8 
65±16 
77±10 

48±8 

42±10 



130±22 

120±35 
110±7 
140±8 
67±15 

82±8 
57±12 

82±15 

220±60 
160±28 
130±38 
130±8 
150±10 
67±18 
92±12 
75±18 



50±18 
130±50 

83±8 
67±17 
130±30 
70±10 

68±5 
180±62 

40±7 
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TABLE 3 

Measured Intensities, Linewidths & Centroids 



GMC 


Transition 


^mol 






Av 






(K km s~^) 


(mK) 


{km s~^) 


(km s~^) 



G 

(02:41:54.12) 
(59:35:59.6) 



H 

(02:41:53.76) 
(59:36:04.7) 



SiO(l-0;v=l) 
SiO(l-0;v=0) 
C2H(|-1) 

HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fe-lfe) 

C34S(2-1) 
HNCO(5o5-4o4) 
SiO(l-0;v=l) 
SiO(l-0;v=0) 
C2H(|-i) 
C2H(i-i) 
HCO+(1-0) 
HNC(l-O) 
HC3N(10-9) 
N2H+(l-0) 
CH30H(2fe-lfc) 

C34S(2-1) 
HNCO(5o5-4o4) 



<19 
1.2±0.5 
6.1±2 
4.1±2 
~15±7 
17±1 
<1.5 
10±1 
25±2 
3.3±2 
33±2 
<19 
<1.0 
<3.3 
<3.3 
<15 
7.8±1 
<1.0 
<1.9 
8.4±2 
<2.2 
;S4.8 



<88 
52±17 
110±17 
51±17 
200±38 
200±16 

<26 
140±15 
300±19 
47±19 
320±29 

<88 

<23 
60±17 

<34 
150±38 
130±16 

<26 
73±15 
200±19 

<28 
110±29 



28±4 
12±4 
27±12 
15±7 
20±2 

23±3 
-27±2 
40±6 
21±2 



5.8±9 

6.2±6 

10±5 

20±6 
-32±2 

11±4 



32±13 
55±10 
72± 28 
62±15 
55±5 

57±7 
55±3 
50±15 
43±3 



78±22 

35±13 
58±17 

57±10 
60±5 

32±10 



Note. — is the main-beam briglitiiess temperature based on a 8'' aperture. Imol is tiie peak integrated intensity in units of K km s^^ for tlie 
resolution given in Table ^ Uncertainties are based on the RMS noise for the temperatures and intensities, and 1(7 from the least-squared gaussian 
fits for the veloeity information. Upper limits represent 2(T values. 




Fig. 6. — Projections of each map onto the first three principal component maps. Left) Projections onto the first and second principal 
components. Right) Projections onto the second and third component, while the figure plane geometry is such that it may be seen as a 
view looking to the left down the PC 1 axis. 
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Fig. 7. — a) The HCN(1-0)/HCO"'"(1— 0) line ratio. The greyscale ranges from to 1.5 with dark being large ratios and contours are 
in steps of 0.333 for the resolution of the HNC(l-O) dataset. HCN{l-0) is from[MTH08. b) The HCN{l-0)/HNC(l-0) line ratio with 
the greyscale ranging from 0.5 to 2.0 and contours are in steps of 0.45 for the resolution of the HNC(l-O) dataset. c) The HCO+{l- 
0)/N2H"'" (1-0) line ratio with the greyscale ranging from 1.0 to 6.0 and contours are in steps of 1.0. d) The HCO"'"(1-0)/HNC{1-0) line 
ratio with the greyscale ranging from 0.5 to 2.0 and contours are in steps of 0.25. 




Fig. 8.— The Position- Velocity diagrams for a) 13CO(1-0) IIMTHOa) , b) CCH(l-0; 3/2-1/2), and c) HNCO(5o5-4o4). The reference 
locations for each map are within ±1 of 02''41"'55'', 59°36 10 and ±2 km s~^ of -20 km s~^. The position angle is taken at 29°. 
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RIGHT ASCENSION (J2000) 



Fig. 9. — C2H contours overlaid on the 2 cm radio continuum 
ranges from 0.3 mjy bm~^ to 2 mjy bm~^ for a beam of 2. 2 X 1. 2. 



3 119941) . Contours are the same as in Fig. [3]and the greyscale 
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Fig. 10. — The shock tracer peak antenna temperature ratios. Ratios are measured over an 8 aperture centered on each cloud. Downward 
(upward) facing arrows are upper (lower) limits. 
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Shocks 



Fig. 11. — Schematic of the overall chemistry of Maffei 2. Left) The schematic adapted from Fig. 10 of lMTOSi to emphasize the similarity 
between the two nuclei's chemical properties. Right) A zoom in of the central molecular ring region with its star formation and molecular 
outflow. 



TABLE 4 

Molecular Abundances in Maffei 2 





N(H2) 


SiO 


CaH 


CaHcor^ 


HCO+ 


HNC 


HC3N 


N2H+ 


C34S 


CH3OH 


HNCO 


A 


2.8(21) 


<6.8(-10) 


<5.4(-8) 




5.4(-9) 


<1.6(-9) 


<2.4(-9) 


<7.9(-10) 


<3.2(-9) 


<1.8(-8) 




B 


9.7(21) 


<2.0(-10) 


<1.5(-8) 




6.1(-9) 


1.4(-9) 


<6.9(-10) 


6.1(-10) 


<9.2(-10) 


1.6(-8) 


<1.6(-8) 


C 


1.9(22) 


<1.0(-10) 


1.8(-8) 


1.8(-8) 


4.2(-9) 


1.9(-9) 


1.0(-9) 


5.8(-10) 


<4.9(-10) 


4.2(-9) 


<2.4(-8) 


D 


2.9(22) 


1.6(-10) 


2.2(-8) 


4.2(-8) 


6.0(-9) 


2.1(-9) 


7.6(-10) 


3.8(-10) 


<3.0(-ll) 


3.0(-9) 


~l.l(-9) 


E 


3.0(22) 


~4.7(-ll) 


2.5(-8) 


4.8(-8) 


4.1(-9) 


2.2(-9) 


4.7(-10) 


2.6(-10) 


<3.0(-10) 


3.7(-9) 


<1.5(-9) 


F 


2.3(22) 


£8.3(-ll) 


9.6(-9) 


9.6(-9) 


3.5(-9) 


2.3(-9) 


1.2(-9) 


7.4(-10) 


8.7(-10) 


1.0(-8) 


8.7(-9) 


G 


2.0(22) 


~1.1(-10) 


1.2(-8) 


2.3(-8) 


7.2(-10) 


1.3(-9) 


<5.0(-10) 


6.0(-10) 


~6.5(-10) 


2.0(-8) 


1.6(-8) 


H 


8.9(21) 


<2.1(-10) 


<1.7(8) 


<2.0(-9) 


<1.6(-9) 


1.3(-9) 


<7.5(-10) 


<2.5(-10) 


<1.0(-9) 


1.5(-8) 


<5.1(-9) 



Note. — Format for entries are a(b)— ax 10^. Each molecule is based on the assumptions of optically thin line emission with Te^; also 10 
K. Upper limits are 2cr. Uncertainties are dominated by systematics and are at least a factor of 3 (sec text for discussion of uncertainties). H2 
column densities arc based on ^^CO(l-O) emission sampled at the resolution of the HNC data with an abundance of [Hg/^^CO] = 7.06x10^ adopted 
HMTHOgft . 

^ Corrected for finite optical depth based of the fine structure line ratios (Table [Tl. 



TABLE 5 
Other Observed Transitions 



Molecule 




Transition 


V 

{GHz) 


T7716 
{mjy /bm) 


GMCs 


HCO 
HCO 
HN"c 

HCCi='CN/HCi='C2N 


lo.i 
lo,i 


-Oo.o; ^; 1-1 
-Oo.o; ^;0-i 

1-0 

10-9 


86.777 
86.806 
87.091 
90.593/90.602 


A A A A 
to to to to 
to to 










89.222(5) 
90.912(5) 


70±15 
10±6 


C 
F 



Note. — Upper limits arc 2(t. 
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TABLE 6 

Selected Intensities Ratios 



J . HCiV(l-O) CCH(3/2-l/2) HCN(l-O) giVC(l-O) CS(2-1) HCN(l-O) 
HiVC(l-O) CC-H"(l/2-l/2) HCO+(l-7)y giV-tac(iroJ C^*S(2-1) CS(2-1) 

A ■■• ■•• •■• >0.79 >1.4 <1.6 

B 2.4±0.4 ■■■ 0.63±0.1 >4.5 >1.4 4.9±1.5 

C 1.6±0.2 2.8±1.5 O.QOibO.l >6.6 >6.0 2.5±0.8 

D 1.3±0.2 1.5±0.3 1.1±0.2 >10 >14 2.1±0.6 

E 1.3±0.2 1.5±0.3 1.3±0.2 >9.4 >13 2.2±0.7 

F 1.1±0.2 2.8±1.8 1.7±0.2 >16 ~4.1 2.4±0.7 

G 0.9±0.1 1.5±0.7 ~2.2 >11 ~3.8 1.1±0.4 

H l.OiO.l ■■■ ■■■ >4.4 >1.0 <2.4 



Note. — The resolution of the measurements are set by the lowest resolution input maps. The HCN(l-O) is fro m|MTH08l and the CS(2-1) from 
JKuno et al.|[2b08l ). In determining the CS/C'^^S ratio it is assumed that CS is spatially completely resolved by the lKuno et al.l (200i) map. If this 
assumption is incorrect true ratios will be slightly lower than those quoted. Uncertainties reflect absolute flux calibration except for the C2H ratios 
which are based on signal-to-noise. Upper limits arc la. 



TABLE 7 
PCA Correlation Matrix 



Maps 




3MM 


C2H(|) 


HCN 


HCO+ 


HNC 


CH3OH 


HNCO 


N2H+ 


HC3N 


C^^S 


SiO 


C2H(i) 


"CO 


1.00 


























3MM 


0.88 


1.00 
























C2H{|) 


0.85 


0.77 


1.00 






















HCN 


0.93 


0.92 


0.84 


1.00 




















HCO+ 


0.89 


0.88 


0.81 


0.98 


1.00 


















HNC 


0.95 


0.87 


0.84 


0.97 


0.94 


1.00 
















CH3OH 


0.68 


0.47 


0.49 


0.43 


0.34 


0.51 


1.00 














HNCO 


0.52 


0.27 


0.32 


0.21 


0.12 


0.33 


0.86 


1.00 












N2H+ 


0.86 


0.71 


0.67 


0.77 


0.72 


0.82 


0.71 


0.58 


1.00 










HC3N 


0.54 


0.53 


0.41 


0.62 


0.60 


0.62 


0.32 


0.16 


0.63 


1.00 








C34S 


0.20 


0.24 


0.17 


0.18 


0.20 


0.19 


0.14 


0.095 


0.20 


0.15 


1.00 






SiO 


0.38 


0.23 


0.34 


0.34 


0.30 


0.38 


0.35 


0.30 


0.35 


0.31 


-0.023 


1.00 




C2H{|) 


0.64 


0.60 


0.70 


0.66 


0.63 


0.63 


0.27 


0.18 


0.45 


0.24 


0.008 


0.32 


1.00 



^ Data from flMTHOa) . 



TABLE 8 
PCA Eigenvectors 



PCA Comp. 


1 


2 


3 


4 


5 


6 


7 


"CO 


0.35 


-0.045 


-0.0043 


0.11 


-0.042 


0.15 


-0.064 


3MM 


0.32 


0.15 


-0.10 


0.12 


-0.086 


0.19 


-0.50 


C2H{|) 


0.31 


0.097 


0.086 


0.24 


0.12 


0.020 


0.81 


HCN 


0.34 


0.21 


0.00 


-0.012 


-0.057 


0.16 


-0.086 


HCO+ 


0.33 


0.28 


-0.029 


-0.024 


-0.055 


0.21 


-0.038 


HNC 


0.35 


0.11 


0.0015 


-0.031 


-0.054 


0.16 


0.038 


CH3OH 


0.23 


-0.55 


-0.03 


0.10 


-0.091 


0.0061 


-0.031 


HNCO 


0.17 


-0.65 


-0.01 


0.17 


-0.050 


-0.10 


-0.061 


N2H+ 


0.32 


-0.19 


-0.089 


-0.12 


-0.20 


-0.061 


0.070 


HC3N 


0.23 


0.067 


-0.13 


-0.68 


-0.29 


-0.52 


0.086 


C34S 


0.080 


-0.0027 


-0.81 


0.0092 


0.56 


-0.11 


-0.0031 


SiO 


0.16 


-0.19 


0.43 


-0.51 


0.66 


0.23 


-0.074 


C2H{i) 


0.24 


0.19 


0.33 


0.37 


0.29 


-0.70 


-0.23 


Egnv. % 


59 


13 


8.3 


6.6 


5.6 


2.7 


1.5 



